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1 INTRODUCTION 



1 Introduction 

IGMtransfer is a numerical code written in Fortran 90/95, intended for simu- 
lating the radiative transfer (RT) of light in the vicinity of the Lya line through 
the intergalactic medium (IGM) . Originally written with the purpose of investi- 
gating how the IGM close to Lya emitting galaxies reshapes the emission line, 
it can also be used to probe the general transmission properties of the IGM, 
thus making it useful in simulations of cosmic reionization. 
Lya is a resonant line, meaning that a photon born in the gas surrounding a hot 
star doesn't escape the galaxy until it has scattered millions of times on neu- 
tral hydrogen, constantly changing direction and frequency. To calculate the 
exact spectrum of Lya photons escaping a galaxy, the full resonant scattering 
RT needs to be solved. Once a Lya photon is the tenuous IGM, the probabil- 
ity of being scattered is much smaller, yet in general not negligible. Although 
the physics of scattering in galaxies and that of scattering in the IGM is not 
inherently different, the difference in physical conditions imposes a natural divi- 
sion of the two schemes: in the dense gas of galaxies, photons are continuously 
scattered in and out of the line of sight, whereas in the IGM, once a photon is 
scattered out of the line of sight, it is "lost", becoming part of the background 
radiation. The probability of a background photon being scattered into the line 
of sight, on the other hand, is vanishingly small. An illustration of this is seen 
in Fig. 1. 

If you haven't already done so, download all the necessary files from the URL 
www.dark-cosmology.dk/~pela/IGMtransfer.html. Besides this documentation, 
the archive IGMtransfer- [version] .tar.gz contains the following files: 



IGMtransf er .f 90 
ProcessIGM.f 90 
F_lam.pro 
fold[IPF] .vim 

test . in 
testdir/ 

- CellData.bin 



- GalData.dat 

tojrmodel . in 
toymodel/ 

- CellData.dat 

- dat2bin.f 90 

- GalData.dat 



Main code. 

Processes the output from IGMtransfer. 
IDL code visuahzing the output from ProcessIGM. 
Three Vim-scripts that structure the contents of the 
above three codes, if your editor is Vim. 
Example input file for IGMtransfer. 
A subdirectory for example files containing the fol- 
lowing two files; 

Example input data file for IGMtransfer, contain- 
ing the physical parameters of the gas in a snapshot 
of a cosmological simulation at z = 3.5. 
Example input file containing physical parameters of 
the galaxies in the snapshot. 
Input file for a small toy model. 
A subdirectory containing the following three files: 
ASCII data file with toy model gas parameters. 
Converts ASCII data to binary, ready for IGM- 
transfer. 

Data for two toy galaxies. 



After the following description of the basic principles and the physics of the 
main code, the individual programs are explained. A more thorough description 
is given in Lamscn ct al. (2011), which represents the work first employing 
IGMtransfer, and in Laursen et al. (2010, my Ph.D. thesis). 



1.1 Underlying concepts 
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Figure 1: Illustration of the difference between the galactic RT and the IGM 
RT. Close to the galaxy, photons are scattered both in and out of the line of sight. 
In the rarefied IGM, photons are mainly scattered out of the line of sight, obviating 
the need for a full Lya RT. The exact value of the distance rg from a galaxy to 
begin the IGM RT is somewhat arbitrary, but is of the order of the virial radius of 
the galaxy. 

1.1 Underlying concepts 

IGMtransfer performs the IGM RT in a "computational box" with a (pos- 
sibly adaptively refined) cell-based structure. The final results are obtained by 
considering the average of the RT performed for many sightlines emerging in 
many directions from many galaxies. For a given sightline, a (normalized) spec- 
trum is emitted, suffering random absorption lines (i.e. the Lya forest) as it is 
continuously redshifted when receding from the galaxy. When the edge of the 
computational volume is reached, the sightline continues in a random, inward 
angle, thus "bouncing" around until the bluest wavelength of the simulated 
spectrum has been redshifted into the Lya resonance (Fig. 2). IGMtransfer 
was originally applied in a non-periodic, spherical volume. The present version 
(vl.l) includes the possibility of using the full volume of the computational box, 
and a future version will include the possibility of utilizing periodic boundary 
conditions, i.e. once the edge is reached, the sightline continues on the opposite 
side of the box. 

The spectrum of each individual sightline is written to an output file which can 
subsequently be processed by the program ProcessIGM. This two-step process 
also enables the user to use the calculated spectra to perform other analyses, 
e.g. investigate them for the relative abundance of different absorption systems. 

1.2 Main output 

The final, main output are the two related quantities: 

1. the transmission function F[X), giving the fraction of light transmitted as 
a function of wavelength (as an average over many sightlines cast through 
the cosmological volume), and 

2. the average transmission T of the IGM in a wavelength interval blueward 
of the Lya line. 
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Figure 2: Illustration of how sightlines are cast through the cosmological volume 
(in the case of a spherical subvolume). To sample sufficiently the full solid angle of 
47r around the individual galaxies, n_los (~10^) sightlines are cast from each galaxy 
(of which four are shown here). Each sightline is started at a distance rg — which 
reflects the distance from the center at which the full scattering RT is no longer 
necessary — from the center of a galaxy and traced until the bluest wavelength of 
the emitted spectrum has been redshifted into resonance. When the edge of the 
spherical volume is reached, the ray "bounces" back, i.e. continues in a random 
inward angle. 

1.2.1 Transmission function 

The resulting value of F{X) at wavelength A for a given sightline is 



The optical depth r is the sum of contributions from all the cells encountered 
along the line of sight: 



Here, riHi.i is the density of neutral hydrogen in the i'th cell, is the distance 
covered in that particular cell, v\\^i is the velocity component of the cell along the 
line of sight, and a{X) is the cross section of neutral hydrogen. Due to the reso- 
nant nature of the transition, the largest contribution at a given wavelength will 
arise from the cells the velocity of which corresponds to shifting the wavelength 
close to resonance. However, for sufRciently high column density absorbers (the 
so-called "damped Lya absorbers) the damping wing of the profile may cause 
absorption at velocities quite far from this. 

If dust is present in the simulation, nHiO'(A) is replaced by nHiCr(A) + ndcrd(A). 
Assuming a spectrum of light Jem (A) escaping a galaxy (simulated, e.g., using 
Monte Carlo simulations, as in Laursen et al. (2009a,b)), the final, observed 
spectrum Jobs (A) is then 



(1) 



cells 




(2) 



Jobs(A) = F(A) X Je,„(A). 



(3) 



1.3 fold [IPF] .vim 
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Figure 3: Illustration of the effect of the IGM on the observed Lyo profile 
emerging from a galaxy at z ^ 3.5. Without taking into account the IGM, the two 
peaks are roughly equally high {left panel). However, when the spectrum is trans- 
mitted through the IGM characterized by the transmission function F{X) {middle 
panel, with the cyangrayish region representing the 68% confidence interval), the 
blue peak is dimished, resulting in an observed spectrum with a higher red peak 
{right panel). The figure is taken from Laursen et al. (2011). 



Figure 3 illustrates this effect. Due to the correlation of the IGM with the 
source, F{X) is non-trivial close to the Lya line, but far from the source, or, 
spectrally speaking, at very blue wavelengths, it becomes a constant function 
of wavelength (as long as one does not enter an appreciably different redshift 
epoch). 

1.2.2 Average transmission 

The transmission function is probably only accurate if you have very high resolu- 
tion in your simulation; if your cell size is of the order of the virial radius of your 
galaxies, you will probably overestimate the absorption. However, IGMtrans- 
FER may still sweeten your life by calculating the average, or rather median, 
transmitted fraction 7" in a large wavelength interval bluward of the Lya line. 
Figure 4 shows such calculated fractions at different redshifts, compared to the 
observations of Songaila (2004). 

1.3 fold [IPF] .Vim 

IGMtransfer is written in one single file. If your editor is Vim, first time you 
open the source code use 

vim -s foldl.vim IGMtransfer . f 90 

This will fold distinct parts of the code into single lines, making it more man- 
ageable. To inspect the lines in a fold, go to that fold and press space. To close 
the fold again, press zc while standing in the fold. If you screw something up, 
exit and open one more time with -s f oldl .vim. 

Similarly, ProcessIGM. f 90 and F_lam.pro can be folded with foldP.vim and 
foldF.vim, respectively. 

If the folds are gone when you open the files next time, try adding these lines 
in your . vimrc: 

au BufWinLeave * mkview 

au BufWinEnter * silent loadview 
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Figure 4: Comparison of observations {black data points) and simulations {col- 
ored data points) of the transmitted flux blueward of the Lyo Une as a function 
of redshift. The figure is taken from Laursen et al. (2011), but the observed data 
points are from Songaila (2004). 



2 IGMTRANSFER 
2.1 Input parameters 

The main input parameters are given in an input file; for the example simula- 
tion, this file is called test . in. The individual input parameters are explained 
below and summarized in Table 1. Additionally, two data files containing the 
physical parameters of the gas and the galaxies, respectively, are required. The 
contents of these are explained in Sec. 2.2 and Sec. 2.3, respectively. 



indir, subdir, CellData, and GalData 

The name of the files containing the gas data and the galaxy data are given in 
the keyword CellData and GalData, respectively. These files are supposed to 
be put in a directory, the name of which is given by the keyword subdir which, 
in turn, is a subdirectory of a mother directory given by indir. 

outdir and loutf ile 

The output of IGMTRANSFER is a file with the name given by the keyword 
loutf ile. It is put in a directory, the name of which is also given by subdir, 
but which is a subdirectory of the directory given by outdir. If you want the 
data files and the output to be in the same directory, simply give the same name 
to indir and outdir. Using two different directories, however, may be useful 
if you wish to save the output on a desk which is backed up. 



2.1 Input parameters 
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Poutf ile 

The output from IGMtransfer is processed by the program ProcessIGM. 
While the average transmission T (Sec. 1.2.2) is written to standard output, the 
transmission function F{X) (Sec. 1.2.1) is written in a text file, the name which 
is given by the keyword Poutf ile. This data can be visualized by means of the 
IDL code F_LAM. 

n_write 

The number of sightlines traced between each time IGMtransfer writes its 
output. 

n_los 

The number of sightlines traced per galaxy. 
SpecRes 

Number of bins into which the spectral range is divided. 
BW 

Two-element vector giving the lower and upper values of the wavelength interval 
in Angstrom to be propagated through the IGM. 

BW_stat 

Two-element vector giving the lower and upper values of the wavelength inter- 
val in Angstrom to be used for calculating T- In principle, the lower value, 
BW_stat(l), could be set equal to BW(1). However, due to the rather broad 
damping wing of Lya, as well as the peculiar motion of the gas elements, absorp- 
tion at a given wavelength is caused by gas occupying a quite large range in real 
space. For this reason, even though the spectrum is traced until BW(1) has been 
redshifted into resonance, the average absorption will begin to decrease some 
Angstrom before. Thus, it is probably best to omit the bluest ten Angstrom or 
so, depending on the redshift^. 

The higher value of BW_stat should be sufficiently far from the Lya line center 
that the state of the IGM — i.e. its density, temperature, ionization state, and 
velocity field — is not correlated with the source. If you want to be, say, 1 
Mpc away from your sources at a redshift of 5, then AA = XqH{z = 5) (1 Mpc) 
/c = 2.3 A, i.e. BW_stat(2) ~ 1213 A. 

If you want to compare to the sample of transmissions by Songaila (2004), use 
the interval [1080,1185] A. 

f _rvir 

As discussed in the introduction, the sightlines should be started sufficiently far 
from the center of the galaxies that scattering into the line of sight is much 
less probable than scattering out of the line of sight. Denoting this distance 
To, if ro is too small, or too large, absorption will tend to be over- or under- 
pin real observations, the reason for not using too blue wavelengths is both not to enter the 
Ly/3 resonance and not to enter an appreciably different redshift epoch; these complications 
are of course not present in IGMtransfer. 
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2 IGMTRANSFER 



estimated, depending on whether the immediate surroundings of the galaxies 
on average cause more or less absorption than the general IGM. Of course rg 
depends on the galaxy in question, in particular its size, as well as the general 
state of the circumgalactic gas. However, as argued in Laurscn ct al. (2011), 
when measuring tq in terms of the virial radii of the galaxies, ro for the different 
galaxies become comparable to each other, and through a convergence study it 
was found that, at least at the redshifts probed in that study (z ^ 2.5 to ^6.5), 
setting ro = l.Srvir gave meaningful results. The keyword f _rvir specifies the 
number of virial radii at which to start the sightlines. 

r_eff 

IGMtransfer was originally constructed to perform the RT in a spherical 
subvolume of a box, with the radius of the sphere given by r_ef f e [0, 1], the 
fraction of the half the side length of the box (i.e. the largest possible sphere 
has r_eff = 1). In order to use the full box, simply set r_eff equal to any 
value larger than a/3, so that the sphere encloses the box fully. 

DustType, f _ion 

If dust is included in the simulations, the string DustType should be set to ei- 
ther 'SMC' or 'LMC', depending on whether the dust is modeled as dust in the 
Small or Large Magellanic Cloud (SMC and LMC), respectively. Motivated by 
observations, the dust density is modeled as being proportional to the density 
of neutral hydrogen, plus a fraction of the ionized hydrogen (as well as to the 
metallicity) . This fraction is given by the keyword f _ion. See Sec. 2.2.6 on how 
to incorporate dust in the simulation. 

z, H_0, Omega_M, and Dmega_L 

The redshift of the simulation, and the present values of the Hubble constant (in 
km Mpc~^), matter, and cosmological constant parameters. These values 
are used to convert the wavelength interval into a physical distance. 
Table 1 summarizes the input parameters. 

2.2 Gas data 

The physical parameters of the intergalactic gas elements are given in one- 
dimensional arrays in binary form in the file indir/subdir/CellData, in a 
nested, space-filling hierarchy. Each array constitutes one record; these records 
are displayed in Tab. 2, along with the first record that specifies a few important 
numbers. The subsequent sections explain in further detail what the numbers 
mean. 

2.2.1 N_cells, D_box, ni, nj, and nk 

The first record consists of the following five numbers: 

N_cells is the total numbers of cells (of all refinement levels). 

D_box is the length of the side of the computational box in proper kpc. That 

is, if your simulation is a box of, say, lOh^^ comoving Mpc, then for h — 0.7 a 

snaphot at z = 3.46 will have D_box = (10^ kpc) / 0.7 / (1-f 3.46) = 3203.075 



2.2 Gas data 
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Input parameters 


Parameter 


Explanation 


indir'^ 


Mother directory for input data subdirectory. 


subdir'^ 


Subdirectory containing the data files. 


CellData'= 


Binary file containing physical parameters of the gas. 


GalData'^ 


Text file containing the galaxy data. 


outdir'^ 


Mother directory for output subdirectory. 


loutf ile'^ 


Output file of IGMtransfer, and input file for Pro- 




CESSlGM. 


Poutf ile'^ 


Output file of ProcessIGM. 


n_write' 


Number of sightlines traced between each output. 


n_los' 


Number of sightlines per galaxy. 


SpecRes' 


Spectral resolution in #bins. 


BW^ 


Two-element vector giving the lower and upper limit 




of the wavelength interval for the RT. 


BW.stat"* 


Two-element vector giving the lower and upper limit 




of the wavelength interval in which to calculate T. 


f _rvir 


Number of virial radii from the center of a galaxy, from 




where to begin a given sightline. 




Radius of sphere in which to perform the RT, in terms 




of half the side length of the box, i.e. G [0, 1] for a 




spherical subvolume, or > a/3 to use the full box. 


DustType'^ 


'SMC or 'LMC, depending on the dust type. 


f_ion'i 


Fraction of ionized hydrogen that contributes to the 




dust density (see Sec. 2.2.6). 




Redshift of snapshot. 




Hubble constant in km Mpc""'^. 


Omega_M^ 


Matter density parameter VIm- 


Omega_L'' 


Cosmological constant density parameter J^a- 



Table 1: Explanation of parameters given in the input file test . in. Super- 
scripts and ' indicate character, real (double), and integer, respectively. 



kpc^. 

ni, nj, and nk is the spatial resolution of the base grid. For a regular, unrefined 
grid, N_cells = nixnjxnk. 

2.2.2 LevelString 

As already mentioned, IGMtransfer includes the possibility of performing the 
RT with "adaptive mesh refinement" (AMR); that is, a grid where a number 
of cells are split up into eight cells which, in turn, may be further split up into 
eight cells, and so on, recursively, for an arbitrary number of times. The criterion 
for splitting a cell will usually be a density threshold, but can in principle be 
anything, e.g. density gradient. 

In the case of an adaptively refined grid, an array containing the refinement 



■^The example simulation described in Sec. 5 is a spherical simulation of lOh ^ comoving 
Mpc at 2 = 3.46, but is actually enclosed in a slightly larger box of D_box = 3400 kpc. 
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2 IGMTRANSFER 



Contents of "CellData" 



Record 



Explanation 



N_cells', D_box°, ni% nj', nk' 



LevelString' 
ii_HIString'' 
TString^ 
V_xString'' 



V_yString'' 
V_zString'' 
ZString'' 
n.HIIString'' 



Total number of cells, total length of com- 
putational box in kpc, base grid resolution 
in y-, and z-direction. 
Refinement levels of cells. Optional. 
Neutral hydrogen number density per cm'^. 
Temperature in Kelvin. 
Gas bulk velocity in a;-direction in km s~^, 
in physical coordinates and with respect to 
the center of the box. 
Gas bulk velocity in y-direction. 
Gas bulk velocity in 2-direction. 
Metallicity in terms of Solar. Optional. 
Ionized hydrogen number density per cm'^. 
Optional. 



Table 2: Superscripts and ° indicate integer, double, and single pre- 
cision real, respectively. All arrays (*String) are of size N_cells. 



level of each cell must also be given. The refinement of a cell is denoted by 
its refinement level £, where C = Q corresponds to the unrefined base grid. 
Thus, if dxQ is the length of a base cell, a cell refined I times has length dx^jl^ . 
Proceeding first in the z-direction, then in the y-direction, and finally in the 
a;-direction, the one-dimensional array LevelString contains the refinement 
levels of each cell. If a cell is refined, this mapping continues in a recursive 
fashion one level below. Figure 5 illustrates the mapping in 2D, while a 3D 
rendering of a toy grid in seen in Fig. 6. In the directory toymodel, an ASCII 
file toymodel.dat containing the relevant cell data for this toy grid is found, 
along with a snippet dat2bin.f90 that converts the data to binary form, ready 
to be read in by IGMtransfer. This structure results in a unique ordering of 
the cells, and thus no x-, y-, and z-positions are necessary. 
The other physical parameters are given subsequently in separate arrays, and 
in the same order (take a look in IGMtransfer's subroutine ReadData to see 
how the records are read). 

In order to use IGMtransfer with AMR, it must be compiled with the fiag 
"DAMR" (see Sec. 2.4). For a regular, non-adaptive grid, omit the LevelString 
record and compile without the DAMR flag. A permitted, but somewhat inef- 
ficient, alternative would be to let LevelString consist solely of nixnjxnk 
zeros. 

2.2.3 n_HIString 

Number density nui of neutral hydrogen in cm~^. 

2.2.4 TString 

Temperature of the gas in Kelvin. Note that in version 1.0 this record contained 
instead the corresponding frequency Dopplcr width. 



2.2 Gas data 
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Figure 5: Two-dimensional example of how the AMR grid is represented in 
the code, with base grid resolution ni = nj = 3, and maximum refinement level 
3. The space-filling curve goes through the grid first in the y-direction, and then 
in the x-direction. The order of cells are given by the indicated numbers. In the 
illustrated case, the array LevelString would consist of the numbers 
[0 0002222111001111 0]. 

2.2.5 V_[xyz] String 

The arrays V_xString, V_yString, and V_zString contain the bulk velocity 
of the gas element in the x-, y-, and z-direction, respectively, in km s"^. The 
bulk velocity is the sum of the elements' peculiar motion and the overall Hubble 
expansion, and are expressed in physical coordinates with respect to the center 
of the box. Note that in version 1.0 this record contained instead the velocities 
in terms of the thermal velocity broadening of the Lya line. 

2.2.6 ZString and n_HIIString 

IGMtransfer includes an optional prescription for dust. The dust cross sec- 
tion is expressed as a cross section (Td(A) per hydrogen nucleus, with the wave- 
length dependence given by a fit to extinction curves of dust in the SMC and the 
LMC by Pci (1992), Wcingartncr & Drainc (2001), and Gncdin ct al. (2008)^ 
Thus, in principle the "dust density" rid should actually be the hydrogen den- 
sity riH, such that the optical depth of dust, as seen by a photon of wave- 
length A traveling a distance r through a cell of uniform gas and dust density, 
is Td — rnn o'd(A). However, since in general the metallicity in a given cell will 
differ from that of the reference metallicity Zq (= ^smc or ^lmc)i i^d is scaled 

■^IGMtransfer uses a faster functional form, seen in the subroutine XsecD, which is ac- 
curate to the sub-%0 scale around the Lya line, and to a few percent in the range [800,2000] 
([800,1500]) A for SMC (LMC) dust. 
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Figure 6: Three-dimensional toy model of a grid with five levels of re- 
finement. Taking a look in the files toymodel.dat and dat2bin.f90 should 
help you figure out the way the data is organized. See also www.dark- 
cosmology.dk/~pela/Phd/AMRmovic.gif for a 3D rotation. 



by the metallicity in that cell, divided by Zq. Moreover, since arguably dust 
tends to be destroyed in regions where hydrogen tends to get ionized, rather 
than scaling with nn, in Laursen et al. (2009b) we argued that a more realistic 
quantity is obtained by scaling with uhi, plus some fraction of rtHn- Thus, the 
"dust density" in the I'th cell is set equal to 



where Zq is Solar metallicity, and < /ion < 1 determines to which degree 
the ionized hydrogen fraction contributes to the dust density: is complete 
destruction, 1 is complete survival; in Laursen et al. (2009b) it is argued that 
setting /ion — 0.01 is a realistic value, and gives Lya escape fractions lying 
roughly midway between using and 1. For more details, see Laursen et al. 
(2009b), where the value of /ion is also thoroughly discussed. 

Thus, to calculate the dust density, two additional records must be provided in 
"CellData", namely the metallicity in terms of Solar (ZString) and the density 
of ionized hydrogen (n_HIIString). Furthermore, the code must be compiled 
with the flag "Ddust" (see Sec. 2.4). Note, however, that since by far most of 
the dust resides in the galaxies, while the largest part of a given sightline goes 
through the IGM, in general the effect of dust is negligible. 



nd,i = 



("-Hl,i + /ion"-Hllj) 



Zj/ZQ 



(4) 



2.3 Galaxy data 
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Contents of "GalData" 
Label Explanation 

X, y, z X-, y-, and z-position of galaxy in kpc, where the 

origin is placed in the middle of the box. Note 
that IGMtransfer transforms coordinates from 
[-D_box/2, +D_box/2] to [0,D_box]. 

r_vir Virial radius in kpc. 

v_x, v_y, v_z a;-, I/-, and z-components of systemic velocity of galaxy 
in km s~^, in the reference frame of the center of the 
box. This velocity will be similar to, but in general 
different from, the bulk velocity of the cell at which 
the galaxy is centered. 

Table 3: Physical properties of the galaxies at which the sightlines are 
started. In the example simulation these were identified as bound structures 
satisfying the following criteria: it must not be a substructure of a larger 
structure, it must have at least 15 star particles, and it must have circu- 
lar velocity larger than 35 km s~^. These criteria were imposed to ensure 
sufficient resolution. 



2.3 Galaxy data 

Every sightline is assumed to initiate just outside a galaxy, where "just outside" 
is specified by the keyword f _rvir, which is the number of virial radii from the 
center. In Laurscn ct al. (2011), we found that the correlation of the velocity 
field and density field with the galaxies causes a non-trivial transmission close 
to the Lya line (seen in Fig. 3), and that starting the sightlines ~1.5 virial 
radii from a galaxy gives meaningful results. The galaxies in those simulations 
each consisted of ^10'^ to ^10'* cells. We emphasize that if your circumgalactic 
environs are not sufficiently well-resolved, you will probably not find the same 
shape of f (A); the average transmission T, on the other hand, is not so sensitive 
to resolution, and does not depend on your chosen value of f _rvir. 
The normalized spectrum is emitted at rest wavelength in the reference frame of 
the center of mass of a galaxy, which in turn may have a peculiar velocity relative 
to the cell at which it is centered. This spectrum is then Lorentz transformed 
between the reference frames of the cells encountered along the line of sight. 
Since the expansion of space is homologous, each cell can be perceived as lying in 
the center of the simulation, and hence this bouncing scheme does not introduce 
any bias, apart from reusing the same volume several times for a given sightline. 
However, since the sightlines "scatter" around stochastically and thus pierce a 
given region from various directions, no periodicities arise in the calculated 
spectra. 

To calculate the above, the appropriate properties of the galaxies are given in 
the text file given by the keyword GalData, consisting of one line for each galaxy 
(plus a "header" of an arbitrary number of lines beginning with "#" ) , and seven 
rows of data, explained in Tab. 3. 

If you are not interested in the properties of the IGM in the vicinity of galaxies 
(F), but merely in the average properties of the IGM (T), a similar file may be 
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3 PROCESSIGM 



constructed, with random initial positions, r_vir = 0, and v_x, v_y, and v_z 
set equal to the distance from the center of the box times the Hubble constant 
H{z). 

2.4 Compilation 

For memory reasons, all parts of the code concerning the AMR structure and 
the dust are included as "preprocessor directives", and are compiled only if 
the appropriate flags arc set. To invoke AMR structuring, compile with the 
flag "-DAMR" (and include a record in the "CellData" flle with the refinement 
levels of the cells); to invoke dust modeling, compile with the flag "-Ddust" 
(and include two records with the metallicity and the ionized hydrogen density, 
respectively) . 

Accordingly, to compile IGMtransfer for performing RT in the test data file 
testdir/CellData.bin, which include both AMR and dust, write 

ifort -03 -fpp -DAMR -Ddust IGMtransf er .f 90 -o IGMtransfer. x 

or 

gfortran -03 -x f 95-cpp-input -DAMR -Ddust IGMtransf er.f 90 -o IGMtransf er.x 

or 

g95 -03 -cpp -DAMR -Ddust IGMtransfer. f 90 -o IGMtransf er.x 

or 

pgf90 -03 -Mpreprocess -DAMR -Ddust IGMtransf er . f 90 -o IGMtransf er.x 

depending on your compiler. The simulation is then run with 

./IGMtreinsf er.x < test. in 

Even if neither of the two optional flags are used, you must still compile with the 
flags that enable the preprocessing; that is, to compile with neither the AMR 

nor the dust option, use 

ifort -03 -fpp IGMtransfer .f 90 -o IGMtransfer .x 
for the ifort compiler, and similarly for other compilers. 

3 ProcessIGM 

The output of IGMtransfer is a huge file containing a spectrum for each 
sightline, i.e. a filesize of 4 x n_gal x n_los x SpecRes bytes. The name of 
this file is given by the keyword loutf ile in the input parameter file. 
To extract F(A) and T{z) from this file, ProcessIGM is called with the same 
input file as was used for IGMtransfer, i.e. 



. /ProcessIGM. X < test. in 
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where the compiled executable is assumed to have been named ProcessIGM.x. 
The transmission function F{X) is determined by calculating in each wavelength 
bin the median of all values of e""^ for all sightlines. To know the dispersion of 
the transmission for different sightlines, also the 16th and the 84th percentiles 
are calculated. The result is written in a formatted file, the name of which 
is given by the keyword Poutfile in the input file, with one row for each 
wavelength bin, and the columns giving wavelength, median, 16th, and 84th 
percentile. The fifth column is not really used, but gives the average in each 
bin. 

The average transmission T is determined by calculating first for each sightline 
j the total transmitted flux in the wavelength interval given by the keyword 
BW_stat, i.e. 

^^g-r(.'thbin)^ (5) 
^^bins 

where ii (12) is the bin number corresponding to the value of BW_stat(l) 
(BW_stat(2)) and nbins = «2 — ii + 1 is the number of bins in the applied 
interval. Subsequently, the median (and 16 and 84 percentiles) of all sightlines 
is calculated. The resulting three numbers are the three last numbers written 
to standard output when running ProcessIGM. 



4 F_LAM 



The output file from ProcessIGM, given by the keyword Poutfile can be 
visualized using the IDL code F_lam, which as input takes the same input file 
as the two previous codes: 

IDL> .r F_lam.pro 

IDL> F_lam, 'test. in' [, outf ile=outf ile] [, xrange=xrange] [, yrange=yrange] 

The three optional arguments are: 

outf ile: A string giving the name of the figure produced by F_lam. De- 
fault value is 'F_lam.eps'. 

xrange: A two-element vector giving the range in Angstrom of the plotted 
wavelength interval. Default value is [1213,1218]. 

yrange: A two-element vector giving the ?/-range. Default value is [0,1.1]. 

If everything has gone well, the output of F_lam should be an . eps file looking 
something like the middle panel of Fig. 3. 



5 Example simulation 

The included example files in the subdirectory testdir is a snapshot at z ^ 
3.5 from (an improved version of) the hydro/gravity-simulations described in 
Sommcr-Larscn ct al. (2003) and Sonmicr-Larscn (2006). This simulation uses 
a TreePM/smoothed particle hydrodynamics-technique rather than AMR. Ac- 
cordingly, the physical parameters of the particles have first been interpolated 
onto an adaptive mesh, using the appropriate smoothing kernels. A Haardt & 
Madau (1996)-like UV background has been assumed, but initiating at z = 10 
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5 EXAMPLE SIMULATION 



-8.5 -6.8 -5.1 -3.4 -1.7 0.0 




Figure 7: Surface brightness map of Lyo emission in the cosmological volume 
of the example simulation, i.e. how the volume would look if the hya radiation 
escaped directly towards the observer. 



rather than z = 6. Furthermore, this simple UV RT scheme has been supple- 
mented by a more elaborate post-processed ionizing UV RT scheme, as described 
in Razoumov & Sommer-Larsen (2006, 2007). 

Figure 7 shows the simulated volume. 



5.1 Toy model 

Additionally, a subdirectory toymodel includes an ASCII file with parameters 
for the toy grid seen in Fig. 6. The purpose of this is primarily to see how the 
data are structured; performing the RT in this is possible, but doesn't really 
give meaningful results. The code dat2bin converts the ASCII data to binary 
form in a file that can be read by IGMtransfer. Note that the data do not 
include dust, so IGMtransfer will have to be compiled with the DAMR flag, 
but without the Ddust flag. 
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7 License and citing 

IGMtransfer and its associated programs are free software, distributed under 
the GNU General Public License, i.e. it may be freely distributed, copied, and 
even modified, as long as any changes in distributed versions are indicated. 
For applications of IGMtransfer leading to publications, please cite Laursen 
et aL (2011). 
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